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Abstract 


We use a recently derived result to extract from two-pion interferometry 
data from pp collisions the radius of the coherent component in the source. 

We find a coherent source radius of about 2/m. 

PACS numbers: 25.75. -q, 25.70. Pq, 25.70 Gh. 

Two-particle Bose-Einstein interferometry (also known as Hanbury Brown-Twiss inten¬ 
sity interferometry) as a method for obtaining information on the space-time geometry and 
dynamics of high energy collisions has recently received intensive theoretical and experimen¬ 
tal attention. Detailed investigations have shown that high-quality two-particle correlation 
data can reveal not only the geometric extension of the particle-emitting source but also its 
dynamical state at particle freeze-out m- 

For a partially coherent source, it was suggested that two-pion interferometry can give 
the coherent degree information of the source 0# In paper jllj, we have generalized 
the two-pion interferometry formula for a partially coherent source and found that two- 
pion interferometry was also senstive to the size of the coherent source component. The 
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physical reason is clear: The particle spectrum distribution is the Fourier transformation of 
the source distribution, the coherent source radius information must hide in the spectrum 
distribution and also must appear in the two-pion interferometry formula. We will shown 
below that a general two-pion interferometry formula satisfy the needs to reveal both the 
degree of coherence and space-time information of the coherent source which was neglected 
in previous publications PJTTHTl . 


For a partially coherent source, one can decompose the pion current into a coherent and 
a chaotic term 
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J(x) = J c {x ) + Jl(x), 


( 1 ) 


with 


Jj(x) = I d 4 x'd 4 pj(x',p)'y(x')exp(ip(x — x')). (2) 

ffere j(x',p ) is the probability amplitude of finding a pion with momentum p , emitted by 
the emitter at x'. y(x') = exp(i(p(x')) is a random phase factor which has been taken away 
from j(x',p). All emitters are uncorrelated when assuming: 


< 7*fab (y) >= ^ 4 fa - v)- 


( 3 ) 


ffere < • • • > means phase average. Then the probability amplitude in momentum space 
reads: 

.7 (p) = J J(x)exp{—ip ■ x)d 4 x 

= J/(p) + J c (p), (4) 


with 


Ji(p) — / Ji(x)exp(—ip ■ x)d 4 x, J c (p) = / J c {x)exp(—ip ■ x)d 4 x. 


Using the relationship 


< Jj(x)J c (y) >= 0, 


( 5 ) 

( 6 ) 
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we have 


<JKp)l(p)>= 0 . ( 7 ) 

Then the single pion spectrum can be expressed as: 

Pi(p) = < J*(P)J(P) > 

= <j*(p)Ji(p)>+j* c (p)j c (p). ( 8 ) 

Using relationship 


< r r*( x ih*(yih( x 2h(y2) > = < 7*{x l )^(x 2 ) >< i*(yi)i(y2) > 

+ < 7*0ub(2/2) >< 7*(yi)l(x2) >, (9) 

we obtain 

< Jl(Pl)Jl(P2)Jl(P2)Jl(Pl) > = < \MPl)\ 2 > • < \Jl{P2 )\ 2 > +| < JKPi)MP2) > I'- (10) 

Performing the phase average and using Eq.(4), Eq.(8) and Eq.(10), we can write the two- 
pion spectrum as 

^2 (Pi 1P2) = < J*(Pl)J*(P2)J(P2)J(Pl) > 

= Pi(pi)Pi(p 2 ) + | < Ji(pi)Jl{p 2 ) > | 2 
+2 Re(< J*M)UP 2 ) > r c (P 2 )Jc(Pi)). ( 11 ) 


Similar to ref. [TO], we dehne the following functions: 


gi(x, k) = J d 4 y < -Jj {x + y/2)J I (x 
g c (x, k ) = J d A yJ*(x + y/2)J c (x 


y/2) > exp (ik ■ y ), 
y/2)exp(ik-y). 


( 12 ) 

( 13 ) 


Here gi(x,p ) and g c (x,p) are Wigner functions which describes the source distribution of the 
chaotic source and the coherent source respectively ||TT)| . [TT)| . [T7)|| . From the above definitions, 
we have: 
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Pi(i,j) = PiiPuPj ) =< Ji(.Pi) J i(Pj ) >= 


gi(x, Pl + Pj ) exp(i(p; - p,) ■ x)d 4 x (14) 


and 


Pc(i,j) = Pc(Pi,Pi) = Jc(Pi)MPi) = 1 9c(x, - t - ) exp(i(Pi - Pj) ■ x)d‘x. 


(15) 


Then the two-pion correlation f un ction for a partially coherent source can be expressed as: 

P-iiPuP-i) 


C 2 (P1,P2) = 


pm ■ PM)’ 


(16) 


with 


Pi(p = -Pi(pi) ■ Pi(pj) + |p/(l> 2)| 2 + 2Re(pi(l, 2)p c (2,1)) 


(17) 


and 


no ) = pi(p,p)+pc(p,p)- 


(18) 


For simplicity [§], we assume that pi(i,j ) and j c (p ) are real, then Eq. (16) can be re-written 


C 2 (pi,P 2 ) = 1 + e (Pi) e (P 2 )Hj(pi,p 2 ) + 2i/e(pi)e(p 2 )(l - e(pi))(l - €{p 2 ))H I {p 1 ,p 2 ) (19) 


with 


Hiipupz) = 


Pi{ 1,2) 


n e(p) = 


pi(p,p) 


( 20 ) 


/ Pl(Pl,Pl)Pl(P2,P2) Fl ^) 

Here e(p) describes the chaotic degree of the source. For a totally chaotic source, e(p) = 1 
while for a totally coherent source, e(p) = 0. Due to the limited statistics, the dependence of 
e(p) on momentum p is often neglected, then the above formula can be simplified as 0. ] 


C 2 (pi,p 2 ) = 1 + e 2 Hj(p 1 ,p 2 ) + 2e(l - e)H I (p 1 ,p 2 ). 


( 21 ) 


Eq. (21) is now widely used in the literature. In the following, we will show that Eq. 
(21) is too simple to approximate Eq.(19) and the information about the radius of the 
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coherent source is lost in Eq.(21). We will derive a formula which is better than Eq.(21) to 
approximate Eq. (19). Eq. (16) can be re-written as 


C 2 {pi,p 2 ) = 1 + 


|P/(1,2)I 2 PKk) 


+ 2 Re 


Vc(l, 2)p/(2,1)\ P?(k) 


( 22 ) 


P?(k) Pi<J>i)Pi(P2) V P?(k) ) Pi(j>i)Pi(p 2 ) 

Here k = (pi + p 2 )/2. Because the single particle distribution is exponential form in pp and 
ee + collisions [Q, we have P 2 (k) / P (pi) / P (p 2 ) ~ 1. Under the assumption that pi(i,j) and 
jc(p) are real, the above formula can then be re-expressed as 


C(phP2) = 1 + t 2 (k)R 2 ,(p 1 .p 2 ) + 2e(fc)(l - e(k))R c ( 2, l)i?/(l,2) 


(23) 


with 


Ri(pi,p 2 ) = Pl ^]^\ ,R c (pi,P2) = Pc ^’ 2 ^ 


(24) 


p I (k,k) , ^'' CXI ' 1,1 '* / p c {k, k )' 

One clearly sees the difference between Eq.(23) and Eq.(21), a extra term R c (p\,p 2 ) appears 
in Eq.(23). Similar to ref. |U|, we take both R I (p 1 ,p 2 ) and R c (p\,p 2 ) as Gaussian functions 
and neglect the dependence of e on momentum, then Eq.(23) can be written down as 


C 2 (pi,p- 2 ) = 1 + e 2 exp(-r 2 j ■ q 2 ) + 2e(l - e)exp(-(r 2 + r 2 c ) • q 2 / 2). 


(25) 


Here ry and r c is the radius of the chaotic source and the coherent source respectively. 
Eq.(25) is a limiting form of the general expressions derived in Ref [jllj] for the case of 
a 1-dimensional Gaussian parametrization. In the above derivation, we did not take into 
account the influence of resonance decay on the two-pion correlations which may distort the 
simple Gaussian function assumed above and may also result in an effective double Gaussian 
two-pion correlation function. In paper [|U|], Kulka et al. used the following formula 

C 2 (pi,p 2 ) = 1 + Ai • exp{-R\ ■ q 2 ) + \ 2 exp(-Rl ■ q 2 ) (26) 

to fit their data and found that Ai = 0.30±o;og, ^2 = 0.45 ± 0.13, R\ = 0.63 ± 0.04/m and 
R -2 = 2.11±0.63 fm with y 2 /NDF = 1.09. In the paper [[□[], Kulka et al found that Eq.(26) 
was not consistent with Eq. (21). But if we take those values into Eq.(25), we will find that 
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rj = 0.63 ±0.04 fm ,r c = 2.01 ±0.67 fm and e = 0.54 ±0.07. It is very interesting to notice 
that the size of the coherent source produced in pp collisions is around 2 fm which is larger 
than the radius of chaotic source. Also from the above comparison we found that Eq.(21) 
is too simple to approximate Eq.(19). 

Conclusions: In this paper, it was shown that the widely used two pion interferometry 
formula (Eq.(21)) was too simple to approximate Eq.(19). The information about the radius 
of coherent source was lost in previous publication. Using Eq.(23) to fit the data, we found 
that the radius of coherent source produced in pp collisions was around 2 fm. 
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